Cosmological evolution of global monopoles 
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We investigate the cosmological evolution of global monopoles in the radiation dominated (RD) 
and matter dominated (MD) universes by numerically solving field equations of scalar fields. It is 
shown that the global monopole network relaxes into the scaling regime, unlike the gauge monopole 
network. The number density of global monopoles is given by n(t) ~ (0.43±0.07)/t 3 during the RD 
era and n(t) ~ (0.25±0.05) /t 3 during the MD era. Thus, we have confirmed that density fluctuations 
produced by global monopoles become scale invariant and are given by 5p ~ 7.2(5.0)o- 2 /i 2 during 
the RD (MD) era, where a is the breaking scale of the symmetry. 



PACS: 98.80.Cq 

Since Kibble pointed out that topological defects may 
have been formed as a consequence of cosmological 
phase transitions 0], many cosmological applications 
have been investigated |^|. In particular, the possibil- 
ity is discussed that topological defects should become 
a source to produce density fluctuations responsible for 
the large-scale structure formation and the cosmic mi- 
crowave background (CMB) anisotropy Recently, the 
BOOMERANG experiment || and the MAXIMA exper- 
iment |I| on observations of the CMB anisotropy revealed 
that our universe is flat, which is consistent with the pre- 
diction of inflation. However, they also found a relatively 
low second acoustic peak. One explanation of such a fea- 
ture is a hybrid model, where primordial fluctuations are 
comprised of adiabatic fluctuations induced by inflation 
and isocurvature ones by topological defects j5). In fact, 
topological defects are compatible with the inflationary 
scenario Q. Moreover, deviations from Gaussianity in 
CMB are reported in |Q. Thus, it is still important to 
clarify the dynamics of topological defects and their im- 
plications on cosmology. 

Among many types of topological defects, particular 
attention has been paid to strings. This is mainly because 
it is confirmed that both the local^J l| and the global 
string networks jl0],[ll| relax into the scaling regime, 
where the number of infinite strings per horizon volume 
is a constant irrespective of cosmic time, so that density 
fluctuations produced by strings become scale invariant. 
Thus, strings have been examined as a possible source 
to produce density fluctuations responsible for the large- 
scale structure formation and the CMB anisotropy [Q. 
Furthermore, global strings have been discussed in the 
context of axion cosmology pO|JT^jr^ |. As the result of 
the breakdown of the global Peccei-Quinn (PQ) U(l) 
symmetry, global strings (axionic strings) are formed. 
Then, axions emitted from axionic strings may signifi- 
cantly contribute to the present energy density so that 



* There are still discussions on the energy loss mechanism of 
local strings, gravitational wave radiation via the loop forma- 
tion, or direct heavy particle emission 



they may play the role of the cold dark matter (CDM). 

On the other hand, monopoles have been less investi- 
gated because the formation of local monopoles, generally 
speaking, leads to conflict with observation. If they are 
diluted by inflation |l4j] , annihilated by the Langackcr-Pi 
mechanism fl5| , or swept away by the domain wall [fl6[ , 
the problem can be evaded but local monopoles leave no 
imprint on cosmological evolution. Unlike strings, how- 
ever, the evolution of monopoles has distinct features de- 
pending on whether they are local or global. In fact, 
the long range force works between global monopoles so 
that their annihilations are much more efficient than the 
local counterpart. Therefore, the global monopole net- 
work may go into the scaling regime, where the number 
of global monopoles per horizon volume is a constant ir- 
respective of cosmic time |l7|-|l9[]. 

Bennett and Rhie made the first numerical simula- 
tions and found the tendency that the number of global 
monopoles per horizon volume is a constant though they 
used the nonlinear a model equation to evolve the scalar 
fields Later, Pen, Spergel, and Turok performed 

similar numerical simulations in both the nonlinear a 
model approximation and the full potential Jl9| . (Sec 
also p0|.) However, due to a lack of computer power, 
they can only run a few realizations so that the number 
of global monopoles per horizon volume still has large 
errors and the scaling property cannot be definitely con- 
firmed. Also, as will be shown later, we must pay at- 
tention to the boundary effects, the grid size effects, and 
the total box size dependence. In fact, the monopole can 
disappear under the periodic boundary condition due to 
the boundary effects. 

Now that we have sufficient computer power, we can 
solve the equation of motion of the scalar fields 4> a with- 
out any approximations and run many realizations in or- 
der to investigate and remove the boundary effects, the 
grid size effects, and the total box size dependence so 
that we can ascertain whether the global monopole net- 
work enters the scaling regime. If so, the number density 
of global monopoles is completely determined. In this 
paper, we report the result of our numerical simulations 
for the evolution of the global monopole network. Also, 
a simple analytic estimate is given. 
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As mentioned above, we directly solve the equation of 
motion for scalar fields in the expanding universe. Let 
us consider the following Lagrangian density for scalar 
fields cj) a (x) (a = 1,2,3), 

= \g^rd»r - V eS [<f) a ,T\, (l) 

where is the flat Robertson- Walker metric and the 
effective potential V e g[(j) a ,T] is given by 



Venire] 



\{<f-V 2 ? + -y-if)- (2) 



Here 4> = yfft^ ', r) = a^l - (T/T c ) 2 , and T c = §Vl5er 
is the critical temperature. For T > T c , the potential 
V c ff[(f> a , T] has a minimum at the origin and the 0(3) 
symmetry is restored. On the other hand, for T < T Cl 
new minima </> = 77 appear and the symmetry is broken. 
In this case the phase transition is of second order. 

Then, the equation of motion for 0° in the expanding 
universe is given by 



4> a (x) + 3H(j) a (x) 



1 



rV 2 fW 



Bit) 2 

+\tt 2 (x)- v 2 }r(x)=0, (3) 



where the dot represents time derivative and R(t) is 
the cosmic scale factor. The Hubble parameter H = 
R(t)/R(t) and the cosmic time t are given by 
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1 = CRD 

2H ~ T 2 



H 2 = —^n.T 4 , t = — = ^f (for RD), 
45m 2 y ' ^ ^ 2 



4tt 3 



H 2 = a {T)—^g*T\ t 
45mp! 



2 _ caid 
3H = T 3 / 2 



(forMD), (4) 



where m p \ = 1.2 x 10 19 GeV is the Planck mass, and g* is 
the total number of degrees of freedom for the relativis- 
tic particles. For the MD case, we have defined a{T) 
[a(T) > 1] as a(T) = p mat (T)/ j0rad (T) = a c (T c /T), 
where p ma t(T) is the contribution to the energy den- 
sity from nonrelativistic particles, p ra ,d(T) the contribu- 
tion from relativistic particles at the temperature T, and 
ol c = Pmat(Ic) / prad(Tc) ■ We also define the dimensionless 
parameter £ as 



iRD 




^ (for RD), 
a 



1/2 



^ (for MD). (5) 



10 and 5 to inves- 



In our simulation, we take Crd.md 
tigate C dependence on the result. 

We start the simulations at the temperature Ti — 2T C , 
which corresponds to = t c /4 (RD) and U = t c /{2^/2) 



(MD). Since the 0(3) symmetry is restored at the initial 
time (Ti > T c ), we adopt as the initial condition the 
thermal equilibrium state with the mass, 



(G) 



which is the inverse curvature of the potential at the ori- 
gin at t = U. 

Hereafter we normalize the scalar field in units of tj , 
t and x in units of U, A is set to be A = 0.25, and the 
scale factor R(t) is normalized as R(l) = 1. 

We perform numerical simulations in seven different 
sets of lattice sizes and spacings for the RD case and the 
MD case (See Tables | and |l[). In all cases, the time 
step is taken as St = 0.01. In the typical case (1), the 
box size is nearly equal to horizon volume (-ff -1 ) 3 and 
the lattice spacing to a typical core size of a monopole 
6x ~ 1.0/(V / Acr) at the final time tj. Furthermore, in 
order to investigate the dependence of C, we arrange case 
(7) with C = 5. We have simulated the system from 
10 [(2), (3), (5) and (6)] or 50 [(1), (4), and (7)] different 
thermal initial conditions. Also, in order to investigate 
the effect of the boundary condition (BC), we adopt the 
periodic BC and the reflective BC [V 2 4> a (x) — on the 
boundary] . 

In order to judge whether the global monopole network 
relaxes into the scaling regime, we give time development 
of £, which is defined as 



£ = n(t)t 3 , 



(7) 



where n(t) is the number density of global monopoles. 
Before counting the number of global monopoles in the 
simulation box, we must identify global monopoles. We 
introduce two identification methods. In one method (I), 
we use a static spherically symmetric solution, which is 
obtained by solving the equation 



d 2 ^ 
dr 2 



2d(f> 
r dr 



dV c 



cir [ 



= 0, 



(8) 



where 4> a (r, 9, if) 



r)x a jr with x 1 — r sin 9 cos tp, x 2 = 
and the winding number n = 1. 
The boundary conditions are given by 



r sin tt sin tp, x 3 — r cos 



0(r) -> r), (r -> 00), 
0(0) - 0. 



(9) 



Since spacetime is discretized in our simulations, a point 
with tf> — corresponding to a monopole core is not nec- 
essarily situated at a lattice point. In the worst case, a 
point with = lies at the center of a cube. We require 
that a lattice is identified with a monopole core if the 
potential energy density there is larger than that corre- 
sponding to the field value of a static spherically sym- 
metric solution at r = ^/SSx p hy S /2 [<5x p hys = R(t)Sx]. 
Moreover, in order to reduce the error, we regard the 
identified lattices which are connected as one monopole 



2 



core. In the other method (II), a cubic box is identified 
to include a monopole if all (f> a — (a = 1, 2, 3) sur- 
faces pass through the cubic box. In this method, we 
also regard the identified boxes which are connected as 
one monopole core. In fact, as shown later, the results 
with these two identification methods coincide very well. 

As an example, time development of S;rd in cases (1) to 
(6) under the periodic BC for the RD case is described 
in Fig. |l|. We find that after some relaxation period, 
£,rd becomes a constant irrespective of time for all cases. 
Though all arc consistent within the standard deviation, 
£rd tends to increase as the box size does. This is be- 
cause monopoles annihilate more often than those in the 
real universe under the periodic BC so that monopoles 
annihilate too much for smaller box sizes. On the other 
hand, under the reflective BC, we also find that £rd 
becomes a constant irrespective of time after some re- 
laxation period. £rd is listed in Table Q. In this case, 
contrary to the case under the periodic BC, tends 
to decrease as the box size increases. This is because a 
monopole suffers a repulsive force from the boundary so 
that monopoles near the boundary annihilate less often 
and £,rd takes larger values in smaller-box simulations 
due to the boundary effect. Therefore, the real num- 
ber of the monopole per horizon volume lies in between 
those under the periodic and the reflective BC. From the 
results of the largest-box simulations [case (6)], we con- 
clude that £rd converges to a constant £rb ~ 0.43±0.07 
irrespective of the boundary conditions. Therefore we 
can conclude that the global monopole network relaxes 
into the scaling regime in the RD universe. We also show 
time development of £r£> in case (7) under the periodic 
BC in Fig. |^. First of all, as easily seen, the results with 
two identification methods agree very well. Next, £rd 
asymptotically becomes a constant, 0.36 ±0.01, which is 
consistent with all the above cases with £ = 10 within 
the standard deviation. Hence we can also conclude that 
C does not change the essential result. 

For the MD case, we also find that after some relax- 
ation period, the number of global monopoles per horizon 
volume becomes a constant irrespective of cosmic time 
under the periodic BC expect for cases (1), (2), and (3), 
in which global monopoles annihilate too much due to 
the boundary effect. The trend of the boundary effect is 
the same as that for the RD case. Therefore, from the 
results of the largest-box simulations [case (6)], we con- 
clude that £md converges to a constant (md — 0.25±0.05 
irrespective of the boundary conditions. Thus, we have 
completely confirmed that the global monopole network 
goes into the scaling regime in both the RD and MD 
universes. 

The above results can be understood by the following 
simple discussion. The evolution of the number density 
of global monopoles n(t) is described by the Boltzmann 



equation^] 



dn(t) 
dt 



-P(t)n(t) - 3H(t)n{t), 
n(t) 3mn(t) 



T(t) 



t 



(10) 



where R(t) oc t m , P(t) is the probability per unit time 
that a monopole annihilates with an antimonopole, and 
T(t) is the period it takes for a pair of monopoles at rest 
with the mean separation l{t) to pair annihilate. The 
mean separation l(t) is given by l{t) = R(t)r s = n(t) -1 / 3 , 
where r s is the mean comoving separation. Since a con- 
stant attractive force works between a pair of monopoles 
irrespective of the separation length, we assume that the 
relative velocity between them reaches unity at once and 
that they do not spiral around each other for a long pe- 
riod of time. Then, the period T(t) is given by the fol- 
lowing relation, 



T+t 



dt 



dr. 



(11) 



where tg is the initial time where a pair of monopoles are 
at rest. Then, the period T{t) reads 



T(t) 



1—771 

n(i)V3' 



(for t<,«T). 



(12) 



Inserting this into the Eq. (|10[), the number density n(t) 
takes the following asymptotic value, 



n(t) 



27(1 -m) 6 
t 3 



oc t 



-3 



(13) 



From the above asymptotic form, we find that urd (t) ~ 
0.42/i 3 and n M D(t) ~ 0.04/i 3 . For the RD case, n(t) 
for the simulation and the analytic estimate agrees ex- 
cellently up to the proportional coefficient £,rd- The dif- 
ference of the proportional coefficient £md between the 
simulation and the analytic estimate may come from the 
fact that cosmic expansion is so rapid in the matter dom- 
inated universe that our assumption that the relative ve- 
locity between them reaches unity at once breaks down 
to some extent. 

We summarize our results. By directly solving equa- 
tions of motion for scalar fields in the expanding uni- 
verse, we have confirmed that the global monopole net- 
work enters the scaling regime in both the RD and MD 
universes. The number density of global monopoles is 
given by n(i) = (0.43 ± 0.07) /t 3 for the RD case and 
n(t) = (0.25 ± 0.05)/t 3 for the MD case. Then, density 
fluctuations induced by global monopoles are given by 
Sp ~ m(t)n(t)/t 3 ~ 7.2(5.0)cr 2 /i 2 for the RD (MD) era, 



*A similar discussion was done in [Elf in the 2 + 1 dimension. 
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where m(t) ~ 47rcr 2 n _1 / 3 (i) (x t is the mass of a global 
monopole. 
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FIG. 1. Time development of ^rd in cases (1) to (6) under 
the periodic BC for the RD case. Symbols (□) represent time 
development of £rd- The vertical lines denote a standard 
deviation over different initial conditions. 



Case(7) (periodic BC) 
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FIG. 2. Time development of £rd in case (7) under the 
periodic BC for the RD case. Filled squares represent time 
development of £rd for identification method (I) using the 
spherical symmetric solution, blank triangles for that (II) us- 
ing the (f> a = surface. The vertical lines denote a standard 
deviation over different initial conditions. 
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TABLE I. Seven different sets of simulations for the RD case. 



Case 


Lattice 


Lattice spacing (Sx) 


c 


Realization 


Box sylc/H 1 


S 


S 




number 


[unit = UR(t)] 






(at final time) 


(periodic BC) 


(reflective BC) 


(1) 


128 3 


V3/W 


10 


50 


1 (at 125) 


0.28 ±0.19 


1.00 ± 035 


(2) 


256 3 


V6/20 


10 


10 


1 (at 250) 


0.31 ±0.20 


0.83 ±0.15 


(3) 


256 3 


V3/20 


10 


10 


1 (at 125) 


0.35 ±0.21 


0.71 ±0.33 


(4) 


128 3 


v/3/5 


10 


50 


2 (at 125) 


0.37 ±0.06 


0.71 ± 0.11 


(5) 


256 3 


Vs/w 


10 


10 


2 (at 125) 


0.36 ±0.07 


0.58 ±0.12 


(6) 


256 3 




10 


10 


4 (at 125) 


0.36 ±0.01 


0.50 ± 0.03 


(7) 


128 3 


v/6/10 


5 


50 


1 (at 250) 


0.36 ±0.17 


0.61 ±0.21 



TABLE II. Seven different sets of simulations for the MD case. 



Case 


Lattice 


Lattice spacing (Sx) 


c 


Realization 


Box size/ ' H 1 


t 


t 




number 


[unit = UR(t)] 






(at final time) 


(periodic BC) 


(reflective BC) 


(1) 


128 3 


3(100) 1/3 /256 


10 


50 


1 (at 100) 


Disappearance 


0.75 ± 0.36 


(2) 


256 3 


3(200) 1/3 /512 


10 


10 


1 (at 200) 


Disappearance 


0.77 ±0.10 


(3) 


256 3 


3(100) 1/3 /512 


10 


10 


1 (at 100) 


Disappearance 


0.82 ± 0.50 


(4) 


128 3 


3(100) 1/3 /128 


10 


50 


2 (at 100) 


0.19 ±0.09 


0.42 ±0.10 


(5) 


256 3 


3(100) 1/3 /256 


10 


10 


2 (at 100) 


0.21 ± 0.09 


0.41 ± 0.05 


(6) 


256 3 


3(100) 1/3 /128 


10 


10 


4 (at 100) 


0.20 ± 0.02 


0.30 ± 0.02 


(7) 


128 3 


3(200) 1/3 /256 


5 


50 


1 (at 200) 


Disappearance 


0.44 ±0.03 



G 



